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Dynamical dark matter (DDM) is an alternative framework for dark-matter physics in which 
the dark-matter candidate is an ensemble of constituent fields with differing masses, lifetimes, and 
cosmological abundances. In this framework, it is the balancing of these quantities against each 
other across the ensemble as a whole which ensures phenomenological viability. In this paper, we 
examine the prospects for the direct detection of a DDM ensemble. In particular, we study the 
' constraints imposed by current limits from direct-detection experiments on the parameter space of 

T -_^ i DDM models, and we assess the prospects for detecting such an ensemble and distinguishing it from 

. traditional dark-matter candidates on the basis of data from the next generation of direct-detection 

■ experiments. For concreteness, we focus primarily on the case in which elastic scattering via spin- 

independent interactions dominates the interaction rate between atomic nuclei and the constituent 
particles of the ensemble. We also briefly discuss the effects of modifying these assumptions. 

'. I. INTRODUCTION 

Dynamical dark matter (DDM) p], Q has recently been advanced as an alternative framework for dark-matter 
physics. In this framework, the usual assumption of dark-matter stability is replaced by a balancing between lifetimes 
and cosmological abundances across a vast ensemble of particles which collectively constitute the dark matter. Within 
O J this framework, the dark-matter candidate is the full ensemble itself — a collective entity which cannot be characterized 
in terms of a single, well-defined mass, lifetime, or set of interaction cross-sections with visible matter. As a result, 
cosmological quantities such as the total relic abundance ^ t ot of the ensemble, its composition, and its equation of 
state are time-dependent (i.e., dynamical) and evolve throughout the history of the universe. Moreover, for this same 
reason, DDM ensembles also give rise to a variety of distinctive experimental signatures which serve to distinguish 
them from traditional dark-matter candidates. A number of phenomenological and cosmological consequences to which 
DDM ensembles can give rise were presented in Refs. [2|,[3[, along with the bounds such effects imply on the parameter 
space of an explicit model within the general DDM framework. DDM ensembles can also give rise to characteristic 
signatures at colliders [4], including distinctive imprints on the kinematic distributions of the Standard-Model (SM) 
particles produced in conjunction with the dark-sector fields. 

In this paper, we examine the prospects for the direct detection of DDM ensembles via their interactions with 
atomic nuclei — a detection strategy [5] which has come to play an increasingly central role in the phenomenology of 
most proposed dark-matter candidates (for reviews, see, e.g., Ref. @). Indeed, conclusive evidence of nuclear recoils 
induced by the scattering of particles in the dark-matter halo would provide the most unambiguous and compelling 
signal — and moreover the only non-gravitational evidence — for particle dark matter to date. Data from the current 
generation of direct-detection experiments have already placed stringent constraints on many models of the dark 
sector, and the detection prospects have been investigated for a variety of traditional dark-matter candidates at next 
generation of such experiments. Studies in the context of particular multi-component models of the dark sector have 
also been performed @, H| . 

Here, we shall demonstrate that DDM ensembles can give rise to distinctive features in the recoil-energy spec- 
tra observed at direct-detection experiments, and that these features can serve to distinguish DDM ensembles from 
traditional dark- matter candidates. These features include resolvable kinks in the recoil-energy spectra, as well as 
characteristic shapes which are difficult to realize within the context of traditional models — particularly under stan- 
dard astrophysical assumptions. As we shall demonstrate, these features should be distinguishable for a broad range 
of DDM scenarios at the next generation of direct-detection experiments. Of course, the potential for differentiation 
within the appropriate limiting regimes of DDM parameter space accords with those obtained in previous studies 
of two-component models [8]. However, as we shall demonstrate, the assumption of a full DDM ensemble as our 
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dark-matter candidate leads to many distinctive features which emerge in significant regions of parameter space and 
which transcend those which arise for models with only a few dark-sector particles. 

This paper is organized as follows. In Sect. El we review the general aspects of dark- matter direct detection. We 
discuss how considerations related to particle physics, nuclear physics, and astrophysics impact both the differential 
and total rate for the inelastic scattering of dark-matter particles with atomic nuclei, and examine the properties of 
the recoil-energy spectra associated with traditional dark-matter candidates. In Sect. Hilt by contrast, we investigate 
how these results are modified when the dark-matter candidate is a DDM ensemble, and we compare the resulting 
recoil-energy spectra to those obtained in traditional dark-matter models. In Sect. IIV[ we derive a set of constraints 
on the parameter space of DDM models from current direct-detection data, and in Sect. El we discuss the prospects 
for obtaining evidence of DDM ensembles at the next generation of direct-detection experiments and for distinguishing 
such ensembles from traditional dark-matter candidates. Finally, in Sect. IVI( we summarize our results and discuss 
possible directions for future study. 



II. DIRECT DETECTION: PRELIMINARIES 



We begin our study by briefly reviewing the situation in which a traditional dark-matter candidate x with a mass 
m x scatters off a collection of atomic nuclei. In general, the differential rate (per unit mass of detector material) for 
the scattering of such a dark-matter candidate off a collection of atomic nuclei can be written in the form [6| 

7 t~) ^(0) JOC 

where Er is the recoil energy of the scattered nucleus in the reference frame of the detector, a^ x is the x- nil cleus 
scattering cross-section at zero momentum transfer, p l ° c is the local energy density of F(Er) is a nuclear form 
factor, ttin is the mass of the scattered nucleus TV, jj J jy x = m x rriN / (m x + tun) is the reduced mass of the x-nucleus 
system, and I x (Er) is the mean inverse speed of x in the dark-matter halo for a given Er. This mean inverse speed, 
which encodes the relevant information about the halo- velocity distribution of x, is given by 

I X (E R ) = / ^Md s v , (2.2) 

where F x (v) denotes the distribution of the detector- frame velocities v of the x m the local dark- matter halo and 
where v = \v\. The lower limit v m [ n on v follows from the condition that only those x with velocities in excess of the 
kinematic threshold for non-relativistic scattering 
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u mi n — \j 9 2 > (2*3) 
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can contribute to the scattering rate. Moreover, the halo- velocity distribution ^(v) itself is truncated at |iT+iT e | < v esc , 
where v e is the velocity of the Earth with respect to the dark-matter halo, and where v esc is the galactic escape velocity. 
Indeed, any dark-matter particle with a speed in excess of v esc in the rest frame of the dark-matter halo would likely 
have escaped from the galaxy long ago. 

One of the primary challenges in interpreting direct-detection data is that substantial uncertainties exist in many 
of the quantities appearing in Eq. (|2.ip . It is therefore necessary for one to make certain assumptions about the 
properties of the dark-matter halo, the nuclear form factor, etc., in order to make concrete predictions regarding 
the detection prospects for any given theory of dark matter. Consequently, in this paper, we adopt a "standard 
benchmark" set of well- motivated assumptions concerning the relevant quantities in Eq. (|2. 1|) . 

The first class of assumptions which define our standard benchmark are those related to particle physics. In 
particular, we take the dark sector to comprise a traditional dark-matter particle x which scatters purely elastically 
off nuclei. Moreover, spin-independent scattering is assumed to dominate the total scattering rate. It follows that 
&n x m ay be written in the form 



cr (0) 

a N X 



[Zf PX + (A-Z)f nx ] 2 , (2.4) 



where f PX and f nx are the respective effective couplings of x to the proton and neutron, Z is the atomic number of 
the nucleus in question, and A is its atomic mass. In addition, the interactions between x and nucleons are taken to 
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be isospin-conserving, in the sense that f px = f nx ; hence Eq. (|2.4j) reduces to 



= • (2-5) 



Similarly, it is also useful to define the spin-independent cross-section per nucleon at zero-momentum transfer: 



r (SI) = 4 ^nX f 2 _ (0) Mn X f9 fi x 



where /i nx is the reduced mass of the x-nucleon system. This quantity has the advantage of being essentially inde- 
pendent of the properties of the target material, and therefore useful for comparing data from different experiments. 

The second class of assumptions which define our standard benchmark are those related to the astrophysics of the 
dark-matter halo. In particular, the local dark-matter density is taken to be p}^ « 0.3 GeV/cm and the velocity 
distribution of particles in the dark-matter halo is taken to be Maxwellian. From the latter assumption, it follows 
that the integral over halo velocities in Eq. (j2.2ft is 0, [To| 



, / ^min + Ve \ , / ^min ~ V e \ - v * / v * . 

erf - erf =e t, « ,c/v ° Vm\n < v e 



2(^esc + Ve ~ Vmin) - v * / v * , . (0 7) 



r / ^esc \ c I ^min i ~v e^' 1 ~e 1 mm/ / - r - ^ 

erf ( - erf = e esc/ V esc - < ^min < ^esc + v e 

vq J \ v J VqJtt 



> v e 



where vq ~ 220 km/s is the local circular velocity and 



erf ( ^ _ ^E. e -vlc/vi 



(2, 



is a coefficient which is independent of both time and m x . By contrast, v e is time-dependent and modulates annually 
due to the revolution of the Earth around the Sun. However, in this paper, we focus primarily on the time-averaged 
scattering rate observed at a given experiment. We therefore approximate the expression in Eq. (j2.7p by replacing 
v e with its annual average (v e ) « 1.05 vo in what follows. Finally, the galactic escape velocity is taken to be v esc « 
540 km/s, in accord with the results obtained from the RAVE survey 11 . 

The third class of assumptions which define our standard benchmark are those related to nuclear physics. These 
assumptions are collectively embodied by the nuclear form factor F(Er). In our standard benchmark, this form factor 
is taken to have the Helm functional form [l2[ 



F{E R ) = 3Jl( V2m ^ M l) e-™» E « s2 , (2.9) 
y/2m N E R Ri 

where J\{x) denotes the spherical Bessel function, s « 0.9 fm is an empirically-determined length scale, and R\ = 
V R 2 — 5s 2 , with R= (1.2 fm) x A 1 / 3 . Note that there exist particular values of Er at which F(Er) vanishes in this 
form-factor model (due to the zeroes of the Bessel function), in the vicinity of which results derived using Eq. (|2.9|) 
are unreliable. However, turns out that all such values of Er will lie well outside the range relevant for our analysis. 

Of course, deviations from this standard benchmark can have a potentially significant impact on recoil-energy spec- 
tra. For example, the effects of unorthodox coupling structures [HI, [l4[, more complicated velocity distributions jlEj . 
different values of the local dark-matter energy density [Tfj, and alternative form- factor models [I?) have all been 
investigated in the literature. In this paper, however, we shall concentrate on the effects that arise when a traditional 
dark-matter candidate is replaced by a DDM ensemble and hold all other aspects of our standard benchmark fixed. 

The standard benchmark described above already leads to characteristic spectra for traditional dark-matter can- 
didates. In Fig. [TJ we display a set of recoil-energy spectra associated with the spin-independent scattering of such 
dark-matter candidates off xenon nuclei (left panel) and germanium nuclei (right panel) for our standard benchmark. 
In each panel, the curves shown correspond to several different values of m x , and each of these curves is normalized 

such that a^ x = 1 pb. 

The shapes of the curves shown in Fig. [T] are determined primarily by two physical effects which suppress the 
differential event rate at large Er. One of these effects stems from the distribution of particle velocities in the dark- 
matter halo. Since F{v) falls off exponentially at velocities above \v + v e \ ~ vq, the recoil-energy spectra likewise 
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FIG. 1: A representative set of recoil-energy spectra obtained for a traditional dark-matter candidate % with a mass m x 
scattering off a xenon target (left panel) and a germanium target (right panel). 



experience a similar suppression above Er ~ 2(vo + v e ) 2 fij^ x /mN- Moreover, these spectra are also truncated at 
Er ~ 2(v esc + v e ) 2 n 2 Nx /rriN as a result of the galactic escape velocity. These effects are particularly important for 
m x < 40 GeV and become increasingly pronounced as m x decreases. Indeed, we will see in Sect. IHII that these effects 
turn out to play a critical role in the direct-detection phenomenology of DDM ensembles, precisely because they are 
sensitive to m x . 

The other effect which plays a significant role in determining the shape of the recoil-energy spectra shown in Fig. [T] 
has its origin in nuclear physics. As is evident from Eq. (j2.9p . the nuclear form factor also suppresses the differential 
event rate for large Er. This suppression is particularly acute for heavier nuclei such as xenon (A & 131), and 
considerably less so for lighter nuclei such as germanium (A « 73). Note, however, that this effect is independent of 
m x and depends only on Er and the properties of the target material. Note also that the dip in the recoil-energy 
spectra displayed in the left panel of Fig. [1] around Er ~ 100 keV corresponds to the first zero of F(Er) in Eq. (|2.9p . 
As noted above, this turns out to lie well outside the range of Er values relevant for this analysis. 

We see, then, that there are qualitatively two distinct regimes which describe the differing behaviors of the resulting 
recoil-energy spectra. In the "low-mass regime," these spectra are steeply falling and highly sensitive to m x . By 
contrast, in the "high-mass regime," these curves fall more slowly and are less sensitive to m x . As we shall see in 
Sect. El this distinction will ultimately play a critical role in our analysis. 



III. DIRECT DETECTION OF DDM ENSEMBLES 



Let us now examine how the results obtained in Sect. [Ill are modified in the case in which the traditional dark- matter 
candidate is replaced by a DDM ensemble, with all of the other defining characteristics of our standard benchmark 
held fixed. As discussed in Refs. |l|, [2|, dynamical dark matter is a new framework for dark- matter physics in which 
the notion of stability is replaced by a delicate balancing between lifetimes and abundances across an ensemble of 
individual dark-matter components. As such, this framework represents the most general possible dark-matter sector 
that can be imagined while still satisfying astrophysical and cosmological bounds. Furthermore, as discussed in 
Ref. [l|, dynamical dark matter arises naturally in certain theories involving extra spacetime dimensions, and also 
in certain limits of string theory. It is therefore important to consider how the direct-detection phenomenology of 
DDM ensembles differs from that of traditional dark-matter candidates. Indeed, such a study can be viewed as 
complementary to the collider analysis performed in Ref. Q. 

From a direct-detection standpoint, the most salient difference between a DDM ensemble and a traditional dark- 
matter candidate is that, by definition, a DDM ensemble comprises a vast number of constituent fields Xji each 
with a mass rrij and local energy density p l ° c . In general, since multiple states are present in the dark sector, both 
elastic processes of the form XjN ~~ ^ XjN an d inelastic scattering processes of the form XjN ~~ ^ XkN ', with j ^ k, 
can contribute to the total Xj- nilc l eon scattering rate. In this paper, in accord with the assumptions underlying 
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our standard picture of dark-matter physics, we focus primarily on the case in which elastic scattering provides the 
dominant contribution to the total scattering rate for each Xj- (This occurs generically, for example, in situations in 
which the mass-splittings between all pairs of constituent particles in the ensemble substantially exceed 100 keV). In 
this case, each Xj also possesses a well-defined effective spin-independent coupling f n j to nucleons, and consequently 
a well-defined spin-independent cross-section per nucleon <jf£P = ^njfnj/ 71 - The total differential event rate at a 
given detector is obtained by summing over the contributions from each Xj, each of which is given by an expression 
analogous to Eq. (|2.ip . Thus, for an arbitrary DDM ensemble, this total differential event rate takes the form 

7 7") A ) JOC 

e£^ 2 (^(^)> (3-d 



dE R ^ 2m j /j J 2 N 



subject to the constraint Y^j p]° c ~ Ptot- Note that the nuclear form factor depends only on Er and not on the 
properties of the constituent particle Xj- By contrast, the integral over halo velocities depends non-trivially on rrij 
through the kinematic threshold velocity 



U) = E R m N 

min " V 2^- * (3 - 2) 

Note that this remains true even in the case we consider here, in which the velocity distributions for all Xj are taken 
to be essentially identical. 

The cosmology of DDM models is principally described by two characteristic quantities [1]. The first of these is 
the collective (present-day) relic abundance Q tot of the full DDM ensemble, which is simply a sum of the individual 
abundances Qj of the Xj- The second quantity is 

V (3-3) 

where = max{Oj}; this helps to characterize the distribution of the Qj across the ensemble, and in particular 
represents the fraction of Q tot collectively provided by all but the most abundant constituent. Thus rj = effectively 
corresponds to the case of a traditional dark- matter candidate, while rj ~ 0(1) indicates that the full ensemble is 
contributing non-trivially to f^tot- 

In general, the local energy densities p^ oc of the Xj — which play a crucial role in direct detection — need not have 
any relation to their cosmological abundances Qj. However, in typical cosmological models, the local energy density 
of any particular Xj is approximately proportional to its cosmological abundance — i.e., p^°VPtot ~ Qj/O to t- 
Furthermore, we assume that Q tot ~ ^cdm, so that the DDM ensemble contributes essentially the entirety of the 
cold-dark-matter relic abundance ^cdm^ 2 ~ 0.1131 ± 0.0034 determined by WMAP [l8|. Under these assumptions, 
the differential event rate in Eq. (|3.ip may be written in the form 

dR _ 2/^ (1 _ > , . . Jj(En) ^ (3.4) 



oIEr TTlJlo 



where mo and f n o respectively denote the mass and effective coupling coefficient of the most abundant state xo in 
the ensemble. 

For concreteness, we examine the direct-detection phenomenology of DDM ensembles in the context of a simplified 
DDM model. In this model xo is identified with the lightest state in the ensemble, and the mass spectrum of the Xj 
takes the form 

rrij — mo + j 6 Am (3-5) 

with Am > and S > 0, so that the Xj are labeled in order of increasing mass. Moreover, in this model Qj and f n j 
are each assumed to exhibit power-law scaling with rrij across the ensemble, so that these quantities may be written 
in the form 

= £7q I — 



m 
mj\ 
mo J 



fnj = fno(^{ , (3.6) 
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where a and /3 are general power- law exponents. Note that scaling relations of this form emerge naturally in many 
realistic DDM scenarios p], 0, H| . Also note that the direct-detection phenomenology of DDM ensembles depends on 
the present-day values of the Qj and not how these values have evolved in the past. Thus the decay widths of the 
Xj , although crucial for the balancing of lifetimes against abundances within the DDM framework |l| , play no role in 
direct detection. 

For the purposes of direct detection, our simplified DDM ensemble is therefore characterized by two groups of 
parameters: those (namely mo, and f n o) which describe the properties of the most abundant state in the ensemble 
and which would also be necessary in any traditional dark-matter model, and those (namely Am and the scaling 
exponents a, /?, and S) which describe how this information extends throughout the entire ensemble. This is therefore 
a very compact yet flexible formalism for exploring the ramifications of having an entire DDM ensemble as our dark- 
matter candidate. However, the WMAP constraint on Q tot fixes one of these parameters (most conveniently Qq). 
Thus, the recoil-energy spectra to which our simplified DDM model gives rise are completely determined by a, /?, <S, 

(SI) 

mo, Am, and f n o (or equivalently cr^ ; ). Note also that the last of these parameters determines the normalization of 
the recoil-energy spectrum (i.e., the total event rate), but has no effect on the shape of that spectrum. 

Given the scaling relations in Eqs. (|3.5p and (|3.6p . it is straightforward to rewrite the expressions for Q tot and 77, 
as well as the differential event rate dR/dEf>, in terms of these parameters. Indeed, in the context of our simplified 
DDM model, the present-day values of Q tot and 77 are given by 



fttot = ^0^] (l+j 



s Am 
mo 



V = 1 



gAmY 
too / 



(3.7) 



From a DDM perspective, our primary interest is in situations in which the number of constituent particles in the 
dark-matter ensemble is taken to be large. For this reason, we restrict our discussion to cases in which the sums in 
Eq. (j3.7p are convergent even in the limit in which j — >> 00. Imposing this requirement restricts the purview of our 
analysis to cases in which the condition aS < — 1 is satisfied. Likewise, the expression in Eq. (|3.4j) for the differential 
event rate reduces to 



dhjR 7rmo ^-^ \ mo / 

3 




FIG. 2: Contours of 77 as a function of the scaling coefficients a and 5, derived under the assumption that Qtot ~ ^cdm- The 
left, center, and right panels show results for Am/mo = {1,0.1, 10 -3 }, respectively. Note that as 5m/mo — > for fixed a and 
5, we see that 77 — ^ 1 and the full ensemble provides a increasingly significant contribution to Qtot- 
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In Fig. [2j we provide a series of contour plots illustrating the dependence of r] on the scaling coefficients a and S in 
our simplified DDM model. The left, center, and right panels of this figure display results for Am /mo = {1, 0.1, 10 -3 }, 
respectively. The white region appearing in each plot is excluded by the condition ad < — 1. The qualitative results 
displayed in this figure accord with basic intuition: r] is maximized for values of a and 6 which come close to saturating 
the constraint a5 < —1, and smaller values of the ratio Am /mo for fixed a and S yield larger values of m However, 
the quantitative results displayed in Fig. [2] are less intuitive and quite significant. In particular, we see that r] ~ 0(1) 
over a broad range of a and S values, even in cases in which Am ~ mo- Within this region of parameter space, the 
full DDM ensemble contributes non-trivially to ^tot- 




E R [keV] E R [keV] E R [keV] 

FIG. 3: Recoil-energy spectra associated with DDM ensembles scattering elastically off of a xenon target. In each of the 
three panels shown, we have set a — —1.5, ft — —1, and 5 = 1, while the left, center, and right panels correspond to the 
choices mo = {10,30,100} GeV, respectively. The different curves displayed in each panel correspond to different values 
of Am, and each of these curves has been normalized so that the total event rate for nuclear recoils in the energy range 
8 keV < Er < 48 keV lies just below the current bound from XENON100 data. Note that the Am — ► oo limit indicated by 
the dashed black curve corresponds to a traditional dark-matter candidate with a mass m x = mo. The dotted black horizontal 
line indicates a reasonable estimate of the recoil-energy spectrum for background events at the next generation of liquid-xenon 
detectors. 
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FIG. 4: Recoil-energy spectra associated with DDM ensembles scattering elastically off of a germanium target. The model 
parameters a, f3, and 5 have been assigned as in Fig.[3l and we have likewise set mo = {10, 30, 100} GeV in the left, center, and 
right panels of the figure, respectively. Each of these curves has been normalized so that the total event rate for nuclear recoils 
in the energy range 10 keV < Er < 100 keV lies just below the current experimental bound. The dotted black horizontal line 
indicates a reasonable estimate of the recoil-energy spectrum for background events at the next generation of germanium-crystal 
detectors. 

We now examine the recoil-energy spectra which arise in the context of our simplified DDM model and identify 
characteristic features in these spectra. A representative set of such spectra is shown for a xenon target in Fig. [3] and 
for a germanium target in Fig. 2J In each of the three panels shown in each figure, we have set a = —1.5, j3 = —1, 
and 5 = 1, while the left, center, and right panels in each figure correspond to mo = {10, 30, 100} GeV, respectively. 
These values have been chosen in order to illustrate the different effects to which DDM ensembles can give rise. The 
different curves displayed in each panel correspond to different values of Am. Each of the curves displayed in Fig. [3] 
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has been normalized such that the total event rate for nuclear recoils in the energy range 8 keV < Er < 48 keV is 
R — 1.0 x 10 -4 kg _1 day _1 . Likewise, each of the curves displayed in Fig. 3] has been normalized so that this same 
total event rate is obtained for nuclear recoils in the energy range 10 keV < Er < 100 keV. Note that this rate is 
consistent with current experimental limits on the total event rate for both target materials; moreover, these recoil- 
energy ranges are chosen to coincide with those typically considered at experiments based on these respective target 
materials. For reference, we also include a curve (the dotted black horizontal line) in Fig. [3] indicating a reasonable 
estimate of the recoil-energy spectrum for background events at the next generation of liquid- xenon detectors (to be 
discussed in more detail in Sect. IV]). We likewise include an analogous curve in Fig. |4] indicating a reasonable estimate 
of the background spectrum at the next generation of germanium-crystal detectors. 

The results displayed in Figs. [3] and 3] demonstrate that the recoil-energy spectra associated with DDM ensembles 
and those associated with traditional dark-matter models differ very little for large mo- This reflects the fact that the 
shape of the contribution to the recoil-energy spectrum from any individual constituent particle Xj 1S n °t particularly 
sensitive to rrij for rrij > 40 GeV. Consequently, for mo > 40 GeV, the contributions from all of the Xj m the ensemble 
manifest roughly the same profile, and the shape of the overall spectrum differs little from that obtained in traditional 
dark-matter models. 

By contrast, the discrepancy between the recoil-energy spectra associated with DDM ensembles and those associated 
with traditional dark-matter models can be quite striking for small mo- In particular, two distinctive features emerge 
which serve to distinguish the recoil-energy spectra associated with DDM ensembles from those associated with 
traditional dark-matter candidates in this regime. The first of these is an apparent "kink" in the spectrum which 
arises for mo < 20 GeV and large Am. Physically, this kink occurs because the contribution from xo to the differential 
event rate dominates at small Er, but falls sharply as Er increases. By contrast, the contribution from each of the 
remaining, heavier Xj falls f ar l ess sharply with recoil energy; hence these contributions collectively dominate at large 
Er. The kink represents the transition point between these two Er regimes. Similar kinks also arise, for example, in 
the recoil-energy spectra of two-component dark- matter models [8[ . 

The second distinctive feature which appears in the recoil-energy spectra displayed in Figs. [3] and 2] emerges in 
cases in which mo and Am are both quite small. In this case, a large number of the Xj are sufficiently light that the 
profiles of their individual contributions to the differential event rate depend quite sensitively on rrij. Moreover, these 
individual contributions cannot be resolved for small Am; rather, they collectively conspire to produce an upturning 
(indeed, an upward concavity) of the recoil-energy spectrum at low Er. The characteristic "S" -shaped or "ogee"- 
shaped curve which results from this upturning is most strikingly manifest in the Am = {1, 10} GeV curves in the left 
panel of each figure. This ogee shape is a distinctive feature of DDM ensembles and is difficult to realize in traditional 
dark-matter models or in multi-component dark- matter models involving only a small number of dark-sector fields. 
As we shall see in Sect. El both the kink and ogee features highlighted here can serve to distinguish DDM ensembles 
at future direct-detection experiments. 

IV. CONSTRAINING DDM ENSEMBLES WITH CURRENT DIRECT-DETECTION DATA 

We begin our analysis of the direct-detection phenomenology of DDM ensembles by assessing how current experi- 
mental data constrain the parameter space of our simplified DDM model. The most stringent limit on spin-independent 
interactions between dark matter and atomic nuclei is currently that established by the XENON 100 experiment [19| 
on the basis of 224.56 live days of observation [20[ with a fiducial mass of 34 kg of liquid xenon. Two events were 
observed within the recoil-energy window 6.6 keV < Er < 30.6 keV which passed all cuts, and 1.0 ± 0.2 background 
events were expected. Under the standard assumptions about the velocity distribution of the particles in the dark- 
matter halo outlined in Sect. HU etc., this result excludes at 90% C.L. any dark-matter candidate — be it a traditional 
candidate or a DDM ensemble — for which the total rate for nuclear recoils with Er within this recoil-energy window 
fails to satisfy the constraint 



In any arbitrary DDM model, the total event rate R for nuclear recoils observed at a given detector is obtained 
by integrating the differential rate in Eq. (j3.4p over the range of Er values which fall within the particular energy 
window Efl in < Er < £^ ax established for that detector. The contribution to this total rate from each Xj is a l so 
scaled by an acceptance factor Aj(En) which depends both on its mass rrij and on recoil energy. In our simplified 
DDM model, we therefore have 



R < 4.91 x 10" 4 kg" 1 day 



-l 



(4.1) 



R = 




dE R F 2 (E R ) J2^i(E R )I J (E R ) ll+f 



3=0 



, s Am 



) 



a+2/3-1 



(4.2) 
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In practice, the dependence of Aj(En) on both Er and mj tends to be slight over the range of Er values typically 
considered in noble-liquid and solid-state detectors. We therefore approximate Aj(Ef>) « 0.5 for all Xj m our analysis 
of the XENON100 constraint, in accord with the acceptance values quoted in Ref. |19j . 
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FIG. 5: Contour plots showing the 90% C.L. limit from XENON100 on the spin-independent cross-section per nucleon 
(in cm -2 ) of the lightest constituent particle xo in our simplified DDM model. The panels appearing in the top, center, and 
bottom rows show results for f3 = {0, —1, —2}, respectively. The panels appearing in the left, center, and right columns show 
the results for S = 0.75, 1,2, respectively. In each panel, we have set a = —1.5. 
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In Fig. O we display a series of contour plots showing the 90% C.L. limit in Eq. (|4.1j) expressed as a bound on cr^ 

7 ( S] 



in our simplified DDM model. Of course, for large Am, the 90% C.L. limit value of a nQ approaches the limit [20|] on 



a^x for a traditional dark- matter candidate with mass m x « mo- However, when Am is small and a larger number of 
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states contribute significantly to the total event rate, the experimental limit can differ substantially from that obtained 
for a traditional dark-matter candidate. Such deviations become particularly pronounced for mo < 10 GeV, in which 
case the majority of nuclear-recoil events initiated by the lightest constituent xo in the ensemble have Er values 
which lie below the detector threshold. In this region, the heavier Xj collectively provide the dominant contribution 
to the total event rate. However, it is evident from Fig. [5] that the heavier Xj can a ls° play an important role in the 
direct-detection phenomenology of DDM models even in the regime in which mo > 10 GeV. 

V. DISTINGUISHING DDM ENSEMBLES AT FUTURE DETECTORS 

We now examine the potential for distinguishing DDM ensembles from traditional dark-matter candidates at future 
direct-detection experiments. Of course, an initial discovery of either a DDM ensemble or a traditional dark-matter 
candidate at a given direct-detection experiment would take the form of an excess in the total number of nuclear- 
recoil events observed above the expected background. Our principal aim is therefore to determine the degree to 
which replacing the traditional dark-matter candidate with a DDM ensemble — keeping all other aspects of our 
standard benchmark unchanged — would result in a discernible deviation in the recoil-energy spectra measured at 
such experiments, once such an excess is observed. For concreteness, we consider the situation in which the total 
scattering rate lies just below the sensitivity of current experiments, so that a sizable number of signal events is 
observed. 

Our procedure for comparing the recoil-energy spectrum associated with a given DDM ensemble to the spectrum 
associated with a traditional dark-matter candidate with mass m x is analogous to that used in Ref. [4] to compare 
invariant-mass distributions at the LHC in the corresponding theories. Similar procedures were also used in Ref. js|. 
In particular, we partition each of the two spectra into bins with widths greater than or equal to the recoil-energy 
resolution AEr at the minimum Er in the bin. We then construct the x 2 statistic 

X>x) = Y. [Xk ~ £ \ {mx)? , (5-D 

k 

where the index k labels the bin, is the expected population of events in bin k in the DDM model, Ek(m x ) is 
the expected population of events in bin k in the traditional dark-matter model, and o\ is the variance in X^ due 
to statistical uncertainties. Since the Xk are distributed according to a multinomial distribution, it follows that 
g\ = Xk(l — Xk/N e ), where N e denotes the total number of signal events observed. 

The proper measure of the distinctiveness of the recoil-energy spectrum associated with a DDM ensemble is not 
the degree to which it differs from that associated with a traditional dark- matter candidate with a particular m x , 
but rather from any such dark-matter candidate. Consequently, we survey over traditional dark- matter candidates x 
with different values of m x with all other assumptions held fixed. Note that the total event rate R — and hence also 

(SI) 

the spin-independent cross-section per nucleon <7n X f° r each value m x — is effectively specified by the signal-event 
count iV e ; thus m x is the only remaining parameter over which we must survey. We then take 

Xmin = min{x 2 (m x )} (5.2) 

as our measure of the distinctiveness of the recoil-energy spectrum associated with a given DDM ensemble. We 
then evaluate a statistical significance of differentiation in each case by comparing Xmin t° a X 2 distribution with 
rib — 1 degrees of freedom. Specifically, this is defined to be the significance to which the p- value obtained from this 
comparison would correspond for a Gaussian distribution. 

For this study, we choose not to limit our attention to any particular experiment, either existing or proposed; 
rather, we investigate the prospects for distinguishing DDM ensembles at a pair of hypothetical detectors, each with 
characteristics representative of a particular class of next-generation direct-detection experiments. The first of these is 
a dual-phase liquid- xenon detector with attributes similar to those projected for XENON1T and future phases of the 
LUX experiment. The other is a germanium-crystal detector with attributes similar to those projected for GEODM 
and the SNOLAB phase of the SuperCDMS experiment. For both experiments, we assume five live years of data 
collection time. 

Each of these hypothetical detectors is characterized by its recoil-energy resolution AEr, recoil-energy window, 
signal acceptance, fiducial mass, and the differential event rate associated with the combined background. For our 
hypothetical next-generation xenon detector, we choose a recoil-energy window 8 keV < Er < 48 keV, a fiducial mass 
of 5000 kg, and a signal acceptance Aj{Er) w 0.5 which is independent of both Er and rnj. We model the energy 
resolution AEr of our detector after that obtained for the combined SI and S2 signals at the XENON100 experiment. 
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This energy resolution is determined from measurements of the detector response for a number of 7-ray calibration 
lines at various energies. The result, expressed in terms of electron-recoil-equivalent energy units keV ee , is [2l| 



0.60 x 



Er 
keV e( 



1/2 



keV e 



(5.3) 



The corresponding energy resolution for nuclear recoils is related to this result by an energy-dependent effective 
quenching factor 



1 



(5.4) 



where C G s = Lnr/L^ is the ratio of the light yield for nuclear recoils to that for electron recoils at zero applied 
electric field, and where S ee = 0.58 and S nr = 0.95 are electric-field-scintillation quenching factors which account for 
the effect of the 530 V/cm applied drift field [22j. The energy resolution AEr for nuclear recoils therefore depends 
on the uncertainties in £ e ff, S ee , and S nr . For Er > 3 keV, the uncertainty in C e ^ is approximately independent of 
Er and given by AC e R « 0.01, while the uncertainty in S ee and S nr is negligible. We therefore find that over the full 
recoil-energy window of our hypothetical detector, AEr is given by 



AE R 



0.36 



(Er\ 

VkeV/ 



/ 0.01 

\£ e fi(ER) 



(ErY 

VkeV/ 



1/2 



keV 



(5.5) 



where all energies are expressed in nuclear-recoil-equivalent units. In our analysis, the width of each bin is set equal 
to the value of AEr at the lowest energy in the bin for this detector. 

We model the differential event rate for the combined background at our hypothetical xenon detector after that 
projected for the combined background at XENON1T. This background rate, after the application of event-selection 
criteria (including both a multiple-scatter veto and an S2/S1 cut), is dominated by electron-recoil events, and in 
particular those from 85 Kr and other impurities within the detector volume. The recoil-energy spectrum associated 
with this background is approximately independent of Er, and for a 85 Kr concentration of 0.5 ppt is given by (23^ 



fdR_\ 
\dE R J BG 



7 x 10~ 9 kg-May^keV" 



(5.6) 



For our hypothetical next-generation germanium detector, we likewise choose the fiducial mass to be 5000 kg. 
Moreover, we choose a recoil-energy acceptance window, energy resolution, and signal acceptance comparable with 
that of the CDMS II experiment. The recoil-energy acceptance window for CDMS II is 10 keV < Er < 100 keV, and 
the energy resolution AEr within this range is given by [24| 



AE h 



"•-(Si) 



1/2 



keV 



(5.7) 



We adopt this same acceptance window and energy resolution for our hypothetical next-generation detector. In order 
to avoid issues related to low statistics for this detector, we adopt a binning scheme coarser than its energy resolution 
would in principal allow. In particular, we set the width of each bin equal to the AEr value of our hypothetical 
xenon-based detector at the smallest Er value in the bin. The signal acceptance for CDMS II varies only slightly, 
from a minimum of A ~ 0.25 at Er values near the endpoints E^ in = 10 keV and i^ ax = 100 keV of the recoil-energy 
acceptance window to a maximum of A ~ 0.32 at Er ~ 20 keV [25|, [26|. We therefore once again approximate the 
acceptance as independent of Er and m j5 and take Aj = 0.3 for all Xj f° r our hypothetical detector. 

In order to obtain a realistic recoil-energy spectrum for the combined background at our hypothetical germanium 
detector, we adopt the following procedure. We model the shape of this spectrum after that observed for the CDMS II 
experiment, which is dominated at Er > 10 keV by the contribution from surface events, and adopt a normalization 
such that the total event rate is Rbg ~ 1.0 x 10 -5 kg -1 day -1 . This is a rate comparable to the background-event 
rate estimates for the SuperCDMS detector at SNOLAB. At recoil energies Er > 10 keV, the background rate at the 
CDMS II detector is dominated by the contribution from surface events. For 10 keV < Er < 25 keV, the recoil-energy 
spectrum is well modeled by [27| 



/ dR \ 
[dE^J 



3.3 x 10" 



x e 



-0.05x(£ 



*/ keV ) kg-May^keV -1 . 



(5i 



12 



From this result, we extrapolate the background spectrum over the full recoil-energy window of our detector. 

We consider the situation in which the total rate for nuclear-recoil events engendered by the DDM ensemble lies just 
below the sensitivity of current experiments, in which case the number of signal events observed at the next generation 
of detectors will be substantial. We examine the DDM differentiation prospects at each of our hypothetical detectors 
independently, in isolation, rather than attempting to correlate the results between the two. For concreteness, we 
adopt a benchmark value of N e = 1000 total signal events at each detector. Note that for our chosen running time of 
five live years, this value of N e is consistent with the XENON100 limits discussed in Sect. II VI throughout the entirety 
of the parameter space of our simplified DDM model which we include in our analysis for both of our detectors. 

In Fig. [6] we show how the projected statistical significance of differentiation obtained with N e = 1000 signal events 
at our hypothetical xenon detector varies as a function of the parameters which characterize our simplified DDM 
model. We find that from among these parameters, the significance is particularly sensitive to the values of mo and 
Am; hence we display our results in (mo, Am) space, with a, /?, and S held fixed in each of the panels shown. The 
panels in the left, center, and right columns of the figure correspond respectively to S = {0.75, 1,2}, and the panels 
in the top, center, and bottom rows correspond to f3 = {—2, —1,0}. In all of the panels shown, we have set a = —1.5. 

The results displayed in Fig. [6] are fundamentally due to the interplay between the individual dR/dER contributions 
from two different classes of Xj within a given ensemble: those with masses rrij < 20 GeV (the "low-mass dark- matter" 
regime) and those with rrij > 20 GeV (the "high-mass dark-matter" regime). As illustrated in Fig. [TJ contributions 
to the recoil-energy spectrum from the Xj m the low-mass regime begin to fall precipitously within or below the 
recoil-energy window for our hypothetical detector. Moreover, the value of Er at which this drop-off occurs is quite 
sensitive to mj for Xj m this regime. By contrast, this suppression effect only becomes manifest for the Xj m the 
high-mass regime at Er values far beyond the recoil-energy-acceptance window of our detector. The spectra fall far 
more gradually with Er for these fields, and their overall shape depends far less sensitively on mj. 

This distinction between these two mass regimes plays a critical role in determining the significance with which 
one can distinguish DDM ensembles from traditional dark- matter candidates at direct-detection experiments. For 
example, it implies that any DDM ensemble in which all of the constituents in the ensemble fall within the high- 
mass regime is generically difficult to distinguish from traditional dark-matter candidates which likewise fall in the 
high-mass regime. Indeed, we see in Fig. [6] that it is quite difficult to distinguish a DDM ensemble in situations 
in which mo > 20 GeV, irrespective of the values of /? and 5. A similar behavior is also manifest in the region of 
parameter space within which mo < 5 GeV and Am > 20 GeV, again regardless of f3 and S. This arises because the 
vast majority of nuclear recoils initiated by any Xj with mj < 5 GeV have Er values which fall below the detector 
threshold E^ in . The contribution to the total recoil-energy spectrum for any Xj with mj in this region is therefore 
essentially invisible. Consequently, in cases in which mo < 5 GeV while Am > 20 GeV, only the contributions from 
the Xj m the high-mass regime are evident and the distinguishing power is once again low. 

By contrast, within other substantial regions of the parameter space of our simplified DDM model, the statistical 
significance of differentiation is quite high. For example, in each panel displayed in Fig. [6j there exists a particular 
range of Am values within which a 5<r significance is obtained for 5 < mo < 20 GeV. Within this region, the kink 
behavior evinced in several of the recoil-energy spectra displayed in Fig. [3] can be distinguished. The range of Am 
values within which this is possible depends primarily on the value of j3. When j3 is small and the coupling to the heavier 
Xj is suppressed, the prospects for distinguishing a DDM ensemble on the basis of this feature becomes significant 
when Am is such that the mass of the next-to-lightest constituent xi li es J us t above the threshold mi ~ 20 GeV of 
the high-mass regime. These prospects remain high until Am reaches the point at which the collective contribution 
to the recoil-energy spectrum from the Xj m the high- mass regime falls below the sensitivity of the detector. As f3 
increases, this contribution remains substantial for larger and larger Am. However, increasing f3 also results in this 
contribution becoming sufficiently large for small Am that it overwhelms the contribution from xo and yields an 
overall spectrum indistinguishable from that of a traditional dark-matter candidate with m x in the high-mass regime. 
The consequences of these two effects are apparent in Fig. [6j which shows how the region of elevated significance 
due to the resolution of a kink in the recoil-energy spectrum shifts from 10 GeV < Am < 50 GeV for f3 = —2 to 
approximately 70 GeV < Am < 800 GeV for (3 = 0. 

Another region of parameter space within which kinks in the recoil-energy spectrum frequently lead to an enhance- 
ment in the significance of differentiation is that within which mo < 5 GeV and 7 GeV < Am < 20 GeV. Indeed, 
such an enhancement is evident in many of the panels in Fig. [6l Within this region, mo is sufficiently light that 
the contribution from xo to the recoil-energy spectrum is hidden beneath the detector threshold, mi lies within the 
low-mass region, and all of the remaining mj with j > 2 lie within the high-mass regime. Thus, within this region, 
Xi plays the same role which xo plays in the region of parameter space discussed above. 

In a number of the panels displayed in Fig. [6] — and especially those in which S < 1 — we obtain a sizeable 
significance of differentiation for our DDM ensemble not merely within this region, but over a substantial region of 
the parameter space within which mo, Am < 20 GeV. Indeed, throughout much of this region, there exist multiple Xj 
with closely-spaced mj in the low-mass regime. When this is the case, the corresponding recoil-energy spectrum for the 
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FIG. 6: Contour plots showing the significance level at which the recoil-energy spectrum associated with a DDM ensemble can 
be distinguished from that associated with any traditional dark-matter candidate which gives rise to the same total event rate 
at a hypothetical direct-detection experiment. This experiment is taken to be a liquid-xenon detector with a fiducial volume 
of 5000 kg and characteristics otherwise similar to those of the proposed XENON IT experiment, as discussed in the text. A 
running time of five live years and an event count of N e — 1000 signal events is assumed. 



DDM ensemble assumes the characteristic ogee shape discussed in Sect. IIIII This ogee shape is a distinctive feature 
of DDM scenarios with small Am, and serves as an effective discriminant between such scenarios and traditional 
dark-matter models. 

As is evident from the results shown in Fig. [6j the significance of differentiation depends on 5 in a somewhat 
complicated manner. Broadly speaking, the significance of differentiation obtained for mo, Am < 20 GeV tends to 
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decrease as 5 decreases, especially for large /?. The primary reason for this is that the density of states in the ensemble 
increases rapidly with rrij when S is small, and thus a greater proportion of Q tot is carried by the Xj i n the high-mass 
regime. Provided that f3 is sufficiently large that a sizeable number of these Xj couple to nucleons with reasonable 
strength, their collective contribution to the differential event rate tends to overwhelm that of the Xj i n the low-mass 
regime. Moreover, even when the low-mass constituents do contribute significantly to the overall rate, the curvature of 
the ogee shape becomes less pronounced — and therefore more difficult to distinguish — as f3 increases. By contrast, 
when 5 is large, the lighter Xj carry a greater proportion of l^tot, and their contributions to the recoil-energy spectrum 
are more readily resolved. However, increasing 5 also increases the mass splittings among the lighter Xj- This has the 
dual effect of pushing a greater and greater number of the Xj m to the high-mass regime and making the individual 
contributions of the remaining constituents in the low-mass region easier to resolve. As a result, the broad regions 
of parameter space throughout which a DDM ensemble could be distinguished on the basis of a characteristic ogee 
feature in the recoil-energy spectrum for small 5 are replaced at large 5 by a set of "islands" in which a kink in the 
spectrum is the distinguishing feature. These effects are already apparent in the right column of Fig. [6l 
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FIG. 7: Contour plots showing the significance level at which the recoil-energy spectrum associated with a DDM ensemble 
can be distinguished from that associated with any traditional dark-matter candidate which gives rise to the same total event 
rate at a hypothetical germanium-crystal detector. The colored regions shown correspond to the same significance intervals 
as in Fig. [3] The detector is taken to have a fiducial volume of 5000 kg and characteristics otherwise similar to those of the 
proposed SNOLAB phase of the SuperCDMS experiment, as discussed in the text. Once again, running time of five live years 
and an event count of N e — 1000 signal events is assumed. 



Let us now compare these results to the results we obtain for our hypothetical germanium-crystal detector. In 
Fig. [7| we show the projected statistical significance of differentiation obtained with N e = 1000 signal events at 
this hypothetical detector. The panels in the left, center, and right columns of the figure correspond respectively 
to S = {0.75, 1,2}, and in each of these three panels we have set a = —1.5 and j3 — —1. It is evident from Fig. 
that while the quantitative results obtained for our two hypothetical detectors differ due to differences in the mass 
of the target nucleus, the observed spectrum of background events, etc., the qualitative results are quite similar. For 
5 GeV < mo < 20 GeV, there exists a range of Am values within which the presence of a discernible kink in the 
recoil-energy spectrum leads to a 5<r significance of differentiation. In addition, a similarly high significance is obtained 
for Am, mo < 20 GeV due either to the similar kink features (at large S) or to the characteristic ogee shape to which 
DDM ensembles can give rise when mass splittings are small (at small S). Note that the particular significance values 
displayed in Figs. [6] and [7] depend on the recoil-energy threshold adopted for our hypothetical detectors and on how 
threshold effects are incorporated into the analysis, especially for small mo- However, varying such assumptions does 
not affect the qualitative results of our analysis. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper, we have investigated the potential for discovering a DDM ensemble and differentiating it from a 
traditional dark-matter candidate at the next generation of dark-matter direct-detection experiments. In particular, 
we have assessed the degree to which these two classes of dark-matter candidates may be distinguished on the basis of 
differences in recoil-energy spectra. We have demonstrated that DDM ensembles give rise to a number of characteristic 
features in such spectra, including observable kinks and distinctive ogee profiles. Moreover, we have demonstrated that 
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under standard assumptions, the identification of such features can serve to distinguish a DDM ensemble from any 
traditional dark-matter candidate at the 5<r significance level at the next generation of direct-detection experiments. 
We have found that the prospects for differentiation are particularly auspicious in cases in which the mass splittings 
between the constituent fields in the DDM ensemble are small and in which the mass of the lightest such field is 
also relatively small. Note that this is also a regime in which a large fraction of the full DDM ensemble contributes 
meaningfully to ^cdm- 

It is also interesting to compare the prospects for distinguishing DDM ensembles at direct-detection experiments 
to the prospects for distinguishing them at the LHC. We have demonstrated here that the former are greatest when 
mo < 20 GeV, Am is small, 0.25 < 5 < 2, and the effective couplings between the Xj an d SM particles decrease 
moderately with rrij. By contrast, it was shown in Ref. Q that the latter are greatest when Am is small, 5 < 1, 
and the effective couplings to SM particles increase with mj. Thus, we see that these two experimental methods 
of distinguishing DDM ensembles are effective in somewhat different regions of parameter space, and are therefore 
complementary. However, we note that there is one region in which evidence for a DDM ensemble may manifest itself 
both at direct-detection experiments and at the LHC. This is the region in which mo < 20 GeV, 0.25 < 5 < 0.75, 
the effective couplings to SM particles are roughly independent of mass, and Am is also either quite small or else 
within the range in which an observable kink arises in the recoil-energy spectrum. The simultaneous observation of 
both collider and direct-detection signatures in this case would provide highly compelling experimental evidence for 
a DDM ensemble. 

Needless to say, there are numerous additional directions potentially relevant for direct detection which we have not 
explored in this paper. For example, we have not considered the prospects for distinguishing DDM ensembles at argon- 
or carbon-based detectors, or instruments involving target materials other than xenon and germanium. Likewise, we 
have not considered the prospects for observing an annual modulation in the signal rate at direct-detection experiments 
— a strategy long employed by the DAM A experiment and more recently by CoGeNT. We have also not considered 
directional detection. More generally, we have not considered modifications of the astrophysical assumptions (such 
as the halo-velocity distributions) or nuclear-form-factor model which define our standard benchmark. Finally, we 
have not endeavored to compare or correlate signals from multiple detectors using different target materials. These 
directions are all ripe for further study [28j. 

In a similar vein, in this paper we have restricted our attention to cases in which elastic processes dominate the 
scattering rate for all particles in the DDM ensemble. However, within the context of the DDM framework, inelastic 
scattering processes [29l-l3l| of the form XjN ~~ ^ XkN where j =^ k also occur, and can contribute significantly to this 
rate when Am < (9(100 keV). This possibility is particularly interesting for a number of reasons. For example, in 
DDM scenarios, the final-state particle in such inelastic scattering events can be a heavier particle in the ensemble, as 
in typical inelastic dark-matter models, but it can also be a lighter particle in the ensemble. In other words, inelastic 
scattering in the DDM framework involves both "upscattering" and "downscattering" processes. This latter possibility 
is a unique feature of DDM scenarios, given that the initial-state particle Xj need not be the lightest particle in the 
dark sector. Moreover, as we have demonstrated above, the range of Am relevant for inelastic scattering is also one 
in which the characteristic features to which DDM ensembles give rise are particularly pronounced. 

Some of the consequences of such inelastic processes are readily apparent. For example, let us consider the case in 
which \Srrijk\ <C m^m^, where 5m jk = rrik — mj- Although the matrix element for inelastic scattering is, to leading 
order, of the same form as for an elastic interaction, the kinematics can be very different. In the limit \5mjk\ <C m^, ra&, 
the maximum recoil energy E^ k and minimum recoil energy EJ k possible in inelastic scattering are given by 

E% = ^(l-^±Jl-2^| (6.1) 

where v is the relative velocity of the initial particles. If 5m ?k > 0, then this upscattering process is similar to that 
typically considered in models of inelastic dark matter [29|, |30| , and its basic effect is to narrow the range of recoil 
energies for which scattering is possible for a fixed dark-matter velocity relative to the Earth. A general result of this 
effect is that heavier components of the DDM ensemble are "brought into range" of a direct-detection experiment, 
which can then resolve the recoil-energy endpoint. For 5mjk < 0, however, the range of possible recoil energies is 
broadened. As a result, low-mass members of the DDM ensemble can produce recoils which lie above the recoil-energy 
threshold E^ in of a particular experiment. Moreover, we note that the matrix element for the process XjN ~~ ^ XkN 
determines the matrix element for the process XkN — >• XjN through crossing symmetry. For a DDM ensemble with 
a fixed distribution of densities, the scattering rates of different components can thus be related to each other. All of 
these possibilities will be discussed further in Ref. (32j . 
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